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Abstract. Let a be an arbitrary angle in a random spherical triangle 
A and a be the side opposite a. (The sphere has radius 1; vertices of A are 
independent and uniform.) If some other side is constrained to be vr/2, then 
E(aa) = 3.05.... If instead some other angle is fixed at 7r/2, then E(aa) = 
2.87.... In our study of the latter scenario, both Apery's constant and Catalan's 
constant emerge. We also review Miles' 1971 proof that E(aa) = 7rV2-2 when 
no constraints are in place. 

For any planar triangle, long sides are opposite large angles and short sides are 
opposite small angles. Quantifying this observation for random triangles with either 

• independent Gaussian vertices in the plane, or 

• independent uniform vertices in a compact convex subset of the plane 

seems analytically intractible. We turn attention therefore to random spherical tri- 
angles with independent uniform vertices on the unit sphere. 

A spherical triangle A is a region enclosed by three great circles on the sphere; a 
great circle is a circle whose center is at the origin. The sides of A are arcs of great 
circles and have length a, 6, c. Each of these is < vr. The angle a opposite side a is 
the dihedral angle between the two planes passing through the origin and determined 
by arcs 6, c. The angles /3, 7 opposite sides 6, c are similarly defined. Each of these 
is < TT too. 

Given a random spherical triangle, the univariate density for a is 

- sin(a), < a < TT 

and 

2 

E(a) = -, E(a2) = y - 2. 
Further, a is uniformly distributed on [0, tt] and 

E(«) = -, E(«') = —. 



^ Copyright © 2010 by Steven R. Finch. All rights reserved. 
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It can be shown that a, b, c are independent random variables; hence E(a b) = 7r^/4 
E(ac). In contrast, the density for (a, /3,7) is [1] 

1 sin(/3) sin(7) sin(a)^ 

47r (1 - (cos(/3) cos(7) - sin(/3) sin(7) cos(a))2)^/2 ' 

As special cases, the conditional density for (/3,7) given that a = 7r/2 is 

1 sin(/3) sin(7) 

2^(1 -cos(/3)2cos(7)2)^/'' 

the conditional density for (0,7) given that f3 = Tr/2 is 

1 sin(7) sin(a)^ 

4(1 — sin(7)2 cos(a)2)^^^ 

and the unconditional density for (/3,7) is 

— cos(7) sin(7) + 7 if /3 
TT + cos(7) sin(7) — 7 if /3 

- cos(/3) sin(/3) + /3 if /3 
TT + cos(/3) sin(/3) - /3 if /3 



27r sin(/3)2 sin(7)^ 



7 > and /3 + 7 < TT, 
7 < and /3 + 7 > TT, 
7 < and /3 + 7 < TT, 
7 > and /3 + 7 > vr. 



These facts will be needed later. 



1. Univariate Densities 
Sides a, b, c are pairwise independent; thus the conditional density for b given c = n/2 
remains unchanged (the sine density on [0,7r]). Angles a, /3, 7 are uncorrelated 
but pairwise dependent. Therefore the case of two angles, plus two other scenarios 
involving opposite side and angle, yield interesting new results. 



1.1. Angle P, for Fixed Angle 7. 

7 = 7i/2 is 



The conditional density for /3 given that 



2sin(/3)2 
1 f 



- cos(/3) sin(/3) + (3 if < /3 < n/2, 
TT + cos(/3) sin(/3) - /3 if 7r/2 < /3 < tt 

-cot(/3) + /3csc(/3)2 ifO</3<7r/2, 
cot(/3) + (vr - /3) csc(/3)2 if 7r/2 < /3 < vr. 



It follows that 



where 



E (3 



TT 

^=2 



71 

2' 



E /3 



TC 

^=2 



^-^C(3) 

2 4^ ^ 



C(3) = E 
fc=i 



is Apery's constant 
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1.2. Side c, for Fixed Angle 7. By the Law of Cosines for Sides: 
cos(c) = cos(a) cos(6) + sin(a) sin(6) cos(7) 

we obtain 

cos(c) = cos(a) cos(fe) 

if 7 = 7r/2. Let u = cos(a), v = cos(6), w = uv, z = arccos(w). Then u, v are 
independent uniform on [—1, 1], that is, with density 



f{u,v) = 
By [21 E] , the density of w is 



1/4 if - 1 < M < 1 and - 1 < w < 1, 
otherwise. 



00 1 

/^(*-f)B* = i/^<'"'4* 



where e{w,t) = 1 if — 1 < w/t < 1, e{w,t) = otherwise. We obtain 

, . 1 



\ia\w\ 



Since < z < n and 



the density of z is 



dz 



dw 



a/1 - w"^ sin(z) 



hiz) 



g(cos(z)) 1 . / M I / M 
1 = — sm{z) In cosu • 



sin(2) 



It follows that the conditional density for side c, given 7 = 7r/2, has a singularity at 
c = 7i/2 and 



E c 



vr \ vr 



7 



E c 



vr 



7 



vr 



-6 + — + IG 



where 



is Catalan's constant [5]. 
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1.3. Angle 7, for Fixed Side c. By the Law of Cosines for Angles: 
— cos(7) = cos(q;) cos(/3) — sin(a) sin(/3) cos(c) 

we obtain 

cos(7) = — cos(a) cos(/3) 

if c = 7r/2. Let u = cos(a), v = cos(/3), w = —uv, z = a,Tccos{w). The Jacobian 
determinant of H- {u,v) is 



sinl^a J 






sin(/3) 



sm 



(a)sin(/3) = Vl^^VT 



because 0<a<7r,0</3<7r. Thus u, v have density 

1 



1 y/l-U^VT^V^ 



1 



1 



271 ^yl-U^^/T^ 27r (1 _ „2^2)3/2 " 

By [21 13] , the density of w is 



oo 



w\ 1 



dt 



tJ\t\ 27^fl_^2^3/2 



where e{w,t) = 1 if — 1 < w/t < 1, e{w,t) = otherwise. We obtain 

, , 1 Inlwl 



TT (1 _ ^i;2)3/2 



and hence the density of z is 



hiz) 



g{cos{z)) 



sin(2:) 



1 In I cosf^;) 



TT sm U 



It follows that the conditional density for angle 7, given c = 7r/2, has a singularity at 
7 = 7r/2 and 

This completes our quick survey of univariate densities, for a fixed side or angle. 

2. BivARiATE Moments 
We evaluate E(aa | b = 7r/2) and E(aa | /3 = 7r/2) here, giving precise numerics for 
the former and exact symbolics for the latter. 
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2.1. (Angle a, Side a), for Fixed Side b. The Law of Cosines for Sides: 
cos(a) = cos(6) cos(c) + sin(6) sin(c) cos(a) 

can be expressed as 

w = uv + Vl — u'^Vl — f ^ cos{6) 

where u = cos(6), v = cos(c), w = cos(a), 6 = a. Then u, v, 6 are independent; m, 
f , w are uniform on [—1, 1] in the unconditional case and 9 is uniform on [0, vr]. Fix 
< 6 < 7r/2 for simplicity, then < m < 1. Solving for v in terms of m, w, 6' we 
obtain two solutions 

uw+\ cos(0)|7(l - u^) {v? - w2 + (1 - n2) cos(^)2) 
(p{u,w,6) = 



ipi^u, w, 6) 



M2 + (l-M2)cOs(e)2 

— I cos(6') I — u^) {u^ — ^2 + (1 — u^) cos{6y) 



m2 + (1 - u^) cos{ey 

assuming 

-w^ + {l- u^) cos{ef > 

and, further, 

{w > —u and 6 < n/2) or {w < —u and 6 > n/2) 

for (fi and 

{w > u and 6 < 7r/2) or {w < u and 6 > 7r/2) 
for 'i/'. Observe that the domains for (f, ip overlap when 

{w > u and 9 < 7c/2) or {w < —u and 9 > tt/2), 

that is, the transformation is one-to-one for {w,9) G [— x [0,7r] and two-to-one 
otherwise. Also, the Jacobian determinant of {v, 9) (-> [w, 9) is 

a/1 - u2f cos(6') 
i;, = M . 

V 1 — 

Let 



^{u,9) = ^n2 + (1 -m2)cos(^)2 
for convenience, then E(aa | b) is equal to [2] 
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2-n J J \\6(u,ip{u,w,e),9)\ \d{u,il]{u,w,9),e)\) 



u 



7r/2 — 

In the event h = 7r/2, we have u = 0, 



cos 



V 1 — 



\6{o,ip{o,w,e),e)\ \6{o,ip{o,w,e),e)\ ip\cos{e)\ (-v^)|cos 

ip I cos 

which becomes 



^^'■^ cos(6»)2 ' 2|w| 



^-W^ + C0S(^)2 I COS(^) I y^-w2 + COS(0)2 

and therefore E(a a | 6 = 7r/2) is equal to 

tt/2 cos(9) tt 

1 /■ f 9 w aTccos(w) , ,^ 1 /■ /■ 6'warccos(w) 

' ' ^ ' -.dwdO + - ^ ' =dw dO. 



^{ I cos(^^)y^-«;2 + cos(0)2 cos(0)^-ti;2 + cos(^)2 

This can be reduced to a single integral: 



TT 

[2- 2Fi(l,i,2,cos(^)2)cos(^) 



^ci^ = 3.0538319164380270202505577... 



involving the following Gauss hypergeometric function: 



4 

TT 



7r/2 7r/2 



- / Jl -xfAnitf dt+(l--\ I ^ = dt 



Despite a connection to complete elliptic integrals [6J, this unfortunately seems to b 
as far as we can go. 
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2.2. (Angle a, Side a), for Fixed Angle j3. The Law of Cosines for Angles: 

— cos(q;) = cos(/9) cos(7) — sm(/9) sm(7) cos(a) 
can be expressed as 

w — —uv + y/l — u^y/l — v"^ cos{6) 

where u — cos(/3), v — cos(7), w = cos(q;), 9 — a. Fix < [3 < 7r/2 for simplicity, 
then < M < 1. Solving for v in terms of u, w, 9 we obtain two solutions (p{—u, w, 9), 
il){—u,w,9) as before. Also, the Jacobian determinant of {v,9) ^ {w,9) is 

VI - cos{9) 

S{u, V, 9) = -u . 

Vl — v 

In the event ^ — 7r/2, we have u — Q and an identical formula for |5(0, </?, ^)|~''^ + 
|(5(0, ■0, ^) follows. The distinction with earlier calculations arises from the density 

1^ 1__ _ 1 si^{9f 

4 (rr^r3;^^y^^^(^j^7rT^ ~ 4(i_(i_^;2)cos(^)2f/' 

for (v, 9). Substituting (/? in place of v, we obtain 

1 sin(^)3 1 sin(^)3 



1 



1 ^ cos(6')^ 1 



cos(e)^ cost 



4 (1 _ ^2^3/2 



and therefore E(q! a | /3 = 7r/2) is equal to 



7r/2 cos(6l) 

1 /■ f 9 w Qxccosiw) sin(^)'' 

' ' :dw d9 



2/ I cos(^)y^^i^2q^^^^(l-w;2)^/' 

dw d9. 



TT 

1 /• r 6waiccos{w) sin(^)3 



\% J(.) cos{9)^-w- + cos{9)Hl - w-^'" 



This can be reduced to a single integral: 

TT 

y ^ tan(^) [cos(^) + sin(^) - 1] = ^ [2 + (1 + ln(2)) tt - 46"] 



4 



= 2.8708787614233542583742340... 



using the fact that 

/ arccosfw) — — , dw 

J (1 _ ^2)3/2 ^_^2 + C0S(^)2 



arccos(t(;jy ^ h y -'^ ^ aw 



sin(6')2 \ V 1 — i 1 — w 



,2 
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and 



arccos w 



arccos w 



'-w^ + cos(6')^ 



-w^ + cos{9y 



cos(e) 



if < < 7r/2, 
if 7r/2 < ^ < TT 



cos(0) 



-^cos(^), 



cos(6l) 



cos(6») 



—w'^ + cos(6')^ 



— + cos(6')^ 
1 — w"^ 



dw 



dw 



if < ^ < 7r/2, 
if 7r/2 < ^ < TT 



:i-sin(0)). 



We have not attempted to extend these formulas for P tt/2. It is intriguing that 
quadrantal triangles (6 = 7r/2) should present an unevaluated integral 3.05... while 
right-angled triangles (/3 = vr/2) give an integral 2.87... expressible in closed-form. 



Unconstrained Scenario 



Miles P proved that 



E((a + /3 + 7-7r)(a + 6 + c)) 



-vr 
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where a + /3 + j ~ it is the area V of the spherical triangle and a + b + c is perimeter 
S. (The notation V, S appears to be traditional.) By preceding correlation results, 



3E(aa) + 6 



TT 

T 



3 E(a a) + 6 E(a 6) - 37r E(a) 



-TT^ - 6 



hence E(aa) = 7r^/2 — 2. It remains to verify Miles' argument. 

Up to now, our random spherical triangles have been built using independent 
uniform vertices. From now on, they will be built using independent uniform great 
circles. By duality, E(U S) = 37r^/2 — 6 under either convention. 

Let k independent uniform great circles be placed on the unit sphere. The number 
of polygonal cells determined is k'^ — k + 2 almost always. Randomly select one of the 
cells (endowed with equal weighting) and denote the density for (V, S) by fk{v,s). 
For example, if = 2, then [7] 



/2(t;,s) = -sin - U(s-27r) if < < 27r 
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and S is the Dirac delta function. No formulas for fk{v,s) are known for A; > 3, 
although when k = 3 marginal densities for V and for S are well-understood [6J. 

Let the cells be labeled randomly by the integers 1, 2, 3, . . ., k"^ — k + 2. It 
is not allowed, for example, to specify that cell 1 cover the north pole and that 
cells 2, 3 be adjacent to it. The labeling must be independent of all features of 
the tessellation. Hence, for the preceding experiment, a cell was selected merely by 
generating a uniform integer j G [1, k"^ — k + 2\. This is the most basic sampling 
technique. 

We wish to examine alternative methods for selecting a cell. Suppose that the 
weighting is proportional to cellular area. Let Cj denote the event that a uniform 
point falls in cell j, where 1 < j < A;^ — A; + 2. If the volume Vj of the cell is 
then the probability of Cj is t'/(47r); unconditionally it is Efc(V^)/(47r). The density 
for (V, S) here is 

gk{v, s) = Pfe { Vj- G [v, V + dv] and Sj G [s, s + ds] \ Cj} 

Pfe {Vj G [v,v + dv] and Sj G [s, s + ds] and Cj} 

_ Pfe {Cj \Vj G [f , f + dv] and Sj G [s, s + ds] } fk{v, s) 

{v/{Att)) fk{v,s) ^ vfk{v,s) 
Efe(y)/(47r) Efe(V^) ' 

v'^fk{v,s) 

Suppose instead that the weighting is proportional to cellular perimeter. A uni- 
form great circle hits 2k cells almost always; we then choose one of these cells at 
random. Let denote the event that a uniform great circle hits cell j and cell j is 
subsequently chosen. If the perimeter Sj of the cell is s, then the probability of Cj is 
(s/(27r))(l/(2/c)); unconditionally it is Efc(5')/(47rfc). The density for (V, S) here is 

hk{v, s) = Pfe ^Vj E [v,v + dv] and Sj G [s, s + ds] \ Cjj 

Pfe ^Vj G [v,v + dv] and Sj G [s, s + ds] and Cjj 



thus 



Pfe{q} 

Pfe |Cj \ Vj G [v,v + dv] and Sj G [s, s -|- ds] | fk{v, s) 
{s/{ATck)) fk{v,s) sfkiv,s) 



Ek{S)/{47rk) Ek{S) 
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thus 

^^/^.(-,3) = ^^. (2) 

Here is an equivalent definition of Cj wliicli is more compatible with that of Cj. 
The intersection of two independent uniform great circles (two diametrically-opposed 
points z and —z) hits two cells almost always; we then choose one of these cells at 
random. The new vertex ±z has four new adjacent cells; upon integrating both sides 
of ([1]), it becomes clear that 

4E«(n - 



Ek{V) 



In the same way, with regard to Cj, the new arc forms the boundary between two 
new adjacent cells; upon integrating both sides of ([2]), it becomes clear that 



Therefore 



and, setting k = 3, 



EkiS) 



Ek-iiV) Ek{S) 



2 E2{V) ^ ' 2 Tx 2 2 

as was to be shown. 

For = 3, the number of cellular vertices is 3 almost always. For k = 4, the 
number iV is 3 with probability 4/7 and 4 with probability 3/7. Recursive equations 
in k for second order moments of V, S, N appear in [H |8] which vastly generalize our 
discussion here. 
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5. Addendum 
M. Larry Glasser reduced the integral 3.05... to an expression 
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where 

F f- - 1 1 - - - \ V r(n + l/2)^r(n + 1)^ x" 

' H2'2' ' '2'2'2'^j ~ 8 ^0 r(n + 3/2)3 • 

He and Jonathan Borwein independently found that 

^ ^2 2 2 2 2 / J 2 



where Li2 is the dilogarithm function. Let agm(a;, y) denote the common hmit of 
sequences {a„} and defined via [11] 

ao = X, bo = y, a„ = ^ ^"""^ , 6„ = A/a„_i6„_i for n > 1. 
David Broadhurst's preferred integral for 3.05... is 

^ sin(0) + 0cos(0) 

7r/2 



agm(l, sin(6')) 



because it permits quick high-precision numerical computation 

[1 
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